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Abstract 


This  study  examines  the  temporal  variation  in 
temperature  at  several  levels  from  near  surface  to 
below  the  thermocline  - a depth  of  approximately 
1900  m.  The  time  series  were  collected  over  a 51 
day  period  from  16  February  through  6 April  1976 
from  an  area  approximately  40  Km  northeast  of  the 
island  of  Eleuthera,  Bahamas.  The  study  also  in- 
cludes aui  examination  of  the  temporal  and  depth 
changes  in  the  density  and  Brunt-Vaisala  frequency 
over  the  period  of  the  record.  These  changes  are 
discussed  in  terms  of  observed  changes  in  the  level 
of  activity  of  the  internal  tides. 
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1.0  Introduction  and  Background  for  the  Data 

The  temperature  data  used  in  the  analyses  for  this 
study  were  acquired  during  February  through  April  1976 
using  the  BEAR  Buoy  System.  (BEAR  - acronymn  for 
Bermuda  Eleuthera  Acoustics  Range) . The  system  com- 
prised a semi-taut  surface  buoyed  mooring  system 
using  a thermistor  cable  to  measure  ocean  temperatures 
to  a depth  of  approximately  1900  m.  The  experimental 
site  was  located  approximately  40  Km  northeast  of  the 
island  of  Eleuthera,  Bahamas  (Fig.  1)  at  25°  45’N, 

76°  17'  W.  The  anchor  position  was  situated  on  the 
abyssal  plain  approximately  20  Km  seaward  of  the  base 
of  the  continental  slope  at  a depth  of  4797  m.  Detailed 
descriptions  can  be  found  in  Kronengold  (1976)  , 
Echternacht  (1976)  , and  Echternacht  (1977) . 

The  scientific  intent  of  the  experiment  was  to 
provide  data  to  study  the  following. 

1.  To  examine  the  amplitude  and  relative  changes 
in  amplitude  of  the  internal  tides  along  the 
Eleuthera  shelf  region. 

2.  To  compare  the  results  of  1.  above  with  open- 
ocean  regions;  in  particular,  the  results  from 
MODE  and  IWEX. 

This  report  presents  the  following : 

1.  An  examination  of  temperature  time  series  at 
varying  depths  from  near  surface  to  below  the 
• thermocline.  The  time  series  represent  a 51 
day  period  from  16  February  through  6 April, 


1976 
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2.  An  examination  of  both  depth  and  temporal 
changes  at  the  high  frequency  cutoff  of  the 
internal  wave  spectrum  - the  Brunt  Vaisala. 

3.  A detailed  discussion  of  changes  in  the 
internal  semi-diurnal  tidal  field  over  this 
period. 

4.  A qualitative  discussion  of  observed  meso- 
scale  features  during  the  period  of  record. 

1. 1 The  Data  Acquisition  System 

The  data  were  acquired  using  the  BEAR  Buoy  system. 

The  system,  shown  schematically  in  Fig.  2,  was  comprised 
of  the  surface  buoy,  the  thermistor  cable,  and  the  nylon 
and  ground  tackle  sections.  The  sensor  subsystem  in- 
cluded the  temperature  sensors  and  the  mooring  sensors. 

The  temperature  sensors  (thermistors)  were  encased 
in  the  thermistor  cable.  The  cable,  Fig.  2,  comprised 
the  upper  portion  of  the  mooring  and  consisted  of  an 
inner  core  of  37  conductors  in  a wet  core  configuration. 
The  individual  thermistor  units  were  potted  and  spliced 
into  the  cable.  The  cable,  in  turn,  was  protected  by  a 
double  lay  torque  balanced  steel  armor  cable.  Table  1 
lists  the  nominal  depth  position  of  the  individual  sensor 
units.  Here  the  nominal  depth  is  defined  as  the  reference 
sensor  depth  including  mooring  line  stretch  under  condi- 
tions of  zero  external  forcing. 
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Fig.  2 BEAR  BUOY  SYSTEM 
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The  mooring  sensors  included  a strain  gauge  tensio- 
meter and  two  inclinometers.  The  tensiometer,  mounted 
in  the  clevis  at  the  base  of  the  buoy,  provided  a measure 
of  the  sum  of  the  external  forces  acting  on  the  mooring 
system.  The  inclinometers  were  mounted  on  the  thermistor 
cable:  one  near  the  surface  and  the  other  at  the  lower 
end  of  the  thermistor  cable  (^1900  m depth)  . The  data 
from  the  mooring  sensors  provided  inputs  to  an  analytical 
model  which  was  used  to  determine  the  mooring  configura- 
tion under  varying  conditions  of  external  forcing. 

For  a more  detailed  description  of  the  data  acquisi- 
tion system  the  reader  is  referred  to  Echternacht  (1976) . 
In  that  study  are  included  the  design  specifications 
for  the  buoy  system  and  the  component  subsystems. 

1.2  Methodology  Used  to  correct  the  Data  Base 

Before  the  temperature  data  could  be  used  in  any 
analysis  scheme  it  was  necessary  to  determine  sensor 
depth  locations  at  each  acquisition  time.  As  discussed 
in  many  previous  works  (eg,  Wunsch  and  Dahlen , 1974) 
deep  instrumentation  located  on  surface  buoyed  mooring 
lines  can  undergo  vertical  excursions  as  large  as  several 
hundred  meters.  The  large  changes  in  sensor  depth  result 
from  low  frequency  horizontal  meanderings  of  the  surface 
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Comparison  of  VDCM  vs  XBT  measurements  showed  that 
the  VDCM  reduced  the  error  between  buoy  acquired  and 
XBT  data  by  as  large  as  a factor  of  three  for  condi- 
tions of  moderate  to  strong  external  forcing. 

In  application  the  VDCM  was  run  for  each  acquisi- 
tion time.  The  model  computed  the  sensor  depth  for 
each  sensor  location.  As  detailed  in  Echternacht  (1977) 
the  temperature  signal  was  corrected  for  vertical  dis- 
placement bias  by  subtracting  that  portion  of  the  signal 
due  only  to  the  change  in  mooring  configuration. 
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2.0  The  Temperature  Data 

The  raw  temperature  and  mooring  sensors  data  were 
acquired  using  an  interval  sampling  technique  in  which 
each  sensor  channel  was  sampled  once  at  the  rate  of  one 
channel  per  second  during  the  turn-on  period.  The  in- 
terval between  sampling  periods  was  12  minutes.  The 
raw  data  were  converted  from  the  acquired  digital  units 
to  engineering  units  using  transfer  functions  developed 
for  each  sensor.  The  transfer  functions  included 
specific  calibration  corrections  (refer  to  Echternacht, 
1976).  The  resulting  data  set  yielded  a times  series  with 
occasional  data  drop-outs.  The  latter  was  the  result  of 
equipment  malfunctions  during  the  acquisition  or  telemetry 
stages.  For  purposes  of  this  study  a 51.5  day  "clean" 
section  was  chosen. 

The  51.5  day  (12  minute  interval)  section  was  then 
low-pass  filtered  yielding  time  series  with  the  new  high 
frequency  cut  off  at  3 hours.  The  temperature  data  were 
then  corrected  for  vertical  displacement  using  the  VDCM 
as  discussed  in  Section  1.2. 

The  3-hourly  temperature  time  series  are  given  in 
Fig.  3.  The  data  cover  the  51.5  day  period  from  1200  LST 
Julian  Day  (JD)  46  through  JD  97,  1976  (15  February  - 
6 April) . The  axes  are  Julian  Days  (time)  on  the  abscissa 
and  temperature  (°C)  on  the  ordinate.  The  same  ordinate 
scaling  was  used  for  all  thermistors.  The  records  for 
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sensor  numbers  2 and  3 are  not  included  - those  sensors 
failed  during  the  system  implant.  Each  temperature 
trace  represents  the  temperature  that  occurred  at  the 
corrected  depth  over  the  period  of  the  record.  The 
sensor  numbers  and  depths  appear  to  the  right  of  each 
temperature  trace.  Here  the  corrected  depth  refers  to 
the  nominal  depth  as  given  in  Table  1.  The  obvious 
features  to  be  noted  are  the  following:  (1)  a large 
mesoscale  feature  covering  approximately  a 24  day 
period  from  the  beginning  of  the  record  to  JD  70, 

(2)  semi-diurnal  oscillations.  The  latter  will  be 
discussed  later  on  in  this  report. 

For  a first  qualitative  attempt  to  classify  the 
mesoscale  feature  the  data  were  compared  with  the  Ring 
Criteria  defined  by  Lai  and  Richardson  (1977) . In  that 
study  they  examined  changes  in  isotherm  depth  location 
through  a cold  core  cyclonic  eddy.  Using  the  15°C 
isotherm  as  an  indicator  it  was  noted  that  in  the  north- 
western Sargasso  Sea  the  15°C  location  is  generally  at  a 
depth  of  650  m outside  an  eddy.  At  the  center  of  the 
eddy  the  following  features  are  noted.  (1)  The  15°C 
isotherm  has  been  displaced  upward  by  approximately  500  m 
and  (2)  packing  of  the  isotherms. 

The  data  given  in  Fig.  3 were  low-pass  filtered 
using  a running  mean  of  4 days  and  the  depth  location  of 
selected  isotherms  plotted  for  the  period  of  the  record. 
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The  isotherm  locations  are  given  in  Fig.  4.  As  seen 
in  the  figure  the  15°C  isotherm  is  displaced  by  ap- 
proximately 150  m - much  less  than  that  observed  by 
Lai  and  Richardson.  Secondly,  there  is  little  ob- 
servable isotherm  packing.  The  data  present  evidence 
of  a mesoscale  event  but  fail  the  ring  criteria  test. 
However,  considering  the  proximity  of  the  buoy  site 
to  the  coast  (a/ 40  km)  it  is  obvious  that  an  eddy  core 
would  lie  well  seaward.  Thus  the  only  conclusion  that 
can  be  drawn  is  that  the  event  might  possibly  have  been 
the  edge  of  an  eddy  or  meander  traversing  the  study 
area.  The  importance  of  the  observation,  however,  is 
that  mesoscales  features  are  relatively  uncommon  in  the 


study  region  (25°  45',  76°  17  W) . In  the  work  by  Lai 


and  Richardson  the  number  of  sighted  anomalies  were 
broken  down  by  regions.  For  the  area  30-40  N,  60-80  W 
there  were  948  sightings  as  compared  to  71  for  the  area 
bounded  by  20-30  N,  65-80  W. 
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3.0  Data  Analyses 

In  terms  of  the  3-hourly  temperature  time  series 
presented  in  the  previous  section,  the  intent  of  this 
report  is  to  provide  a detailed  description  of  and  the 
changes  in  the  internal  semi-diurnal  tidal  field  over 
the  period  shown.  From  an  examination  of  the  data  given 
in  Fig.  3,  the  following  are  noted. 

1)  Pronounced  fluctuations  in  temperature  are 
evident  in  the  time  series;  fluctuations  of  a 
semi-diurnal  periodicity.  The  semi-diurnal  signal 
appears  to  be  the  most  energetic. 

2)  The  amplitude  of  the  signal  increases  with  depth 


to  a maximum  at  approximately  1000  m then  decreases 
with  depth  through  and  below  the  thermocline. 

3)  Within  any  given  depth  zone  there  exist  changes 
in  intensity  of  the  signal  with  time  over  the 
51.5  day  period. 

The  distribution  of  energy  in  the  internal  wave  field 

t 

is  a function  of  the  static  stability  of  the  water  column. 
Hence,  the  analyses  to  follow  examine  the  spatial  and 
temporal  changes  in  the  fields  of  temperature,  density, 
sound  velocity,  and  static  stability. 


3.1  Temperature  and  Salinit: 


I 
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one-day  overlap  between  data  sections.  The  filter  was 
applied  to  the  temperature  data  only.  As  discussed  by 
Pickard  (1964)  for  warm  waters  of  low  latitude  origin 
changes  in  the  density  field  are  dominated  by  tempera- 
ture not  by  salinity.  Thus  for  this  study  the  salinity 
is  considered  a conservative  property.  This  assumption 
is  reasonable  considering  the  length  of  the  record. 

The  3-day  mean  temperature  profiles  and  the  mean 
salinity  profile  are  given  in  Fig.  5.  The  mean  salinity 
profile  was  obtained  using  data  from  the  Fleet  Numerical 
Weather  Center  (FNWC) . The  FNWC  data  are  based  on 
climatological  XBT  data  for  the  study  region  during 
the  winter  months.  The  temperature  and  salinity  data 
are  listed  in  Tables  2 and  3 in  the  Appendix. 

The  temperature  profiles,  in  general,  are  typical  of 
the  winter  months  for  the  study  area.  In  comparison  to 
the  summer  regime  (Mooers,  1975)  the  shallow  warm  mixed 
layer  is  absent.  However,  at  the  beginning  of  the  record 
(profiles  1-6)  and  during  the  latter  half  of  the  period 
(profiles  14-25)  there  exists  a thin  layer  of  decreased 
gradient.  For  profiles  1-6  the  change  in  gradient  occurs 
over  the  depth  interval  of  approximately  250  to  550  m and 
over  the  interval  of  250  to  650  m for  the  latter  profiles. 

In  reference  to  Fig.  4 (page  12)  these  correspond  to  the 
periods  JD  47  through  59  and  JD  74  through  97,  respectively. 
The  anomalies  are  due  to  the  mesoscale  feature  seen  in 
Fig.  4. 
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The  salinity  profile  is  typical  for  those  areas 
influenced  by  the  Antilles  Current.  The  characterizing 
feature  is  the  subsurface  salinity  maximum  in  the  region 
of  100  to  300  m depth. 

3.2  The  Density  Field 

The  density  anomaly,  (T^,  was  computed  using  the  3-day 
mean  temperature  data,  the  mean  salinity  profile,  and  the 
corresponding  mean  sensor  depth.  The  as  used  in  this 
study  is  the  in  situ  density  anomaly. 


<Tfe  = ($s  t p - I ) y io  3 O') 

where,  the  subscript  s,t,p  indicates  in  situ. 

(Tt  profiles  were  computed  for  each  mean  data  section. 

Fig.  5 (bottom)  gives  the  profiles  over  the  period  of  the 
record.  The  data  are  listed  in  Table  4 in  the  Appendix. 

In  the  case  of  the  profiles  the  change  in  temperature 

gradient  (Fig. 5,  upper)  due  to  the  mesoscale  feature  results 
in  a noticeable  anomaly  in  0"j..  As  will  be  seen  later  in 
the  discussion  of  the  static  stability  these  data  suggest 
the  possibility  of  an  intrusion  of  a water  mass  of  different 
characteristics . 

It  should  be  noted  that  the  abscissa  scale  for  Figs.  5 
and  6 is  a sliding  scale.  The  vertical  line  separating 
profiles  is  the  minimum  abscissa  value  for  each  curve.  The 
scale  is  given  once  for  each  figure  in  the  upper  left. 
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salinity,  and  depth) . In  the  software  used  the  pressure 
is  obtained  using  the  mean  data  listed  above.  Fig.  6 
(bottom)  gives  the  mean  sound  velocity  profiles  for  each 
mean  section  over  the  51-day  period.  The  data  are  listed 
in  Table  5 in  the  Appendix. 

The  influence  of  the  mesoscale  event  is  quite  evident 
in  the  sound  speed  profiles.  The  depth  level  and  the 
amplitude  of  the  sound  speed  minimum  change  throughout 
the  period.  At  the  beginning  of  the  record  the  level  of 
the  minimum  decreases  in  depth  through  approximately 
profile  6.  During  the  later  half  of  the  record  the  level 
increases  in  depth.  This  strongly  suggests  first  a 
shallowing  followed  by  a deepening  of  thermocline  depth 
during  the  passage  of  the  mesoscale  event.  This  is 
supported  by  the  results  given  in  Fig.  4.  The  amplitude 
of  the  sound  speed  minimum  increases  throughout  the 
period  from  a value  of  1484  ms-1  at  the  beginning  of  the 
period  to  approximately  1488  ms-1  at  the  end  after  the 
passage  of  the  mesoscale  event. 


3. 4 Static  Stability 

Typically  the  ocean  is  characterized  by  stable  density 
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stratification.  Generally,  the  stratification  is  so 
pronounced  that  seawater  at  a given  depth  is  so  much  denser 
than  waters  at  a shallower  depth  that  it  remains  denser 
even  if  raised  to  the  shallower  depth.  The  parameter  used 
to  characterize  the  density  stratification  is  the  Brunt- 
Vaisala  frequency,  N given  by  the  following. 

1,2 <2>  - -£al  - (§-)2  <2> 

where, 

C is  the  local  sound  speed  and 

dft  the  in  situ  density  gradient. 

dz 

2 

In  this  study  N (z)  was  computed  for  each  sensor  depth 
location  using  the  computed  mean  values  of  C and  . The 
density  gradient  was  computed  using  a finite  difference 
scheme . 

2 

An  examination  of  the  N (z)  profiles  over  the  period  of 
record  yields  the  following  (Fig.  6,  upper). 

1.  During  the  beginning  of  the  period  (profiles  1-6)  there 
exist  two  distinct  stable  zones;  the  surface  to  200m  and 
250  to  450m. 

2.  The  level  of  neutral  stability  changes  through  the  period. 
During  the  early  phase  (profiles  1-6)  the  level  is  at  a 
depth  of  approximately  1000m  then  shallows  to  900m 
through  profile  9.  From  profile  9 through  the  remainder 
of  the  period  the  level  gradually  deepens  to  a final 
depth  of  1150m. 
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3.  The  latter  profiles  are  more  typical  of  waters  influenced 
by  the  Antilles  Current;  ie,  the  stable  zone  at  100-300m 
depth. 

These  data  suggest  the  presence  of  two  distinct  water 
masses  during  the  beginning  of  the  period  associated  with 
the  passage  of  the  mesoscale  feature.  This  is  supported 
by  the  perturbation  in  the  level  of  neutral  stability;  ie, 
thermocline  depth.  During  the  earlier  phases  the  fact 
that  the  deeper  stable  zone  maintained  its  identity  and  that 
the  character  of  the  upper  zone  did  not  change  appreciably 
through  the  record  would  suggest  that  the  mesoscale  feature 
did  not  have  a surface  signature  of  any  consequence. 
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4 . 0 Discussion  and  Conclusions 

To  summarize  the  analyses  presented  in  this  study 
the  following  can  be  concluded. 

1.  During  quiescent  periods  (i.e.,  no  mesoscale 
transient  features  present)  the  temperature, 
density,  sound  speed,  and  static  stability 
profiles  are  typical  for  the  study  area  during 
the  winter  period.  The  profiles  exhibit  the 
influence  of  the  Antilles  Current. 

2.  Over  the  period  of  record  fluctuations  of  a semi- 
diurnal periodicity  are  the  dominant  feature  for 
time  scales  less  than  the  local  inertial.  The 
intensity  of  the  semi-diurnal  signal,  moreover, 
was  observed  to  vary  greatly  over  the  51-day 
period.  In  terms  of  the  mode  structure  associated 
with  the  semi-diurnal  wave  field  the  implication 
of  the  above  is  that  the  distribution  of  energy 

by  mode  is  variable  rather  than  constant  as  is 
generally  assumed  for  this  region. 

3.  The  passage  of  the  transient  mesoscale  feature 
through  the  study  area  perturbed,  noticeably,  the 
fields  of  temperature,  density t sound  speed,  and 
static  stability  - particularly  at  intermediate 
depths  between  300  m and  the  thermocline.  The  mode 
shapes  (vertical  displacements  due  to  the  internal 
wave  field)  are  directly  proportional  to  the  Brunt- 
Vaisala  as  given  by  the  following. 
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where , 

Zn  is  the  vertical  displacement, 

O'  the  semi-diurnal  tidal  frequency, 

■P  the  local  inertial,  and 

d.  the  wave  number. 

2 

Obviously,  the  solution  to  Eq.  (3)  using  the  N 
profiles  presented  in  this  study  would  yield  a much  more 
complicated  picture  of  the  internal  wave  field  for  the 
study  area.  In  particular,  the  modal  structure  and  dis- 
tribution of  energy  by  mode  changes  considerably  with  time. 
The  importance  of  these  findings  to  acoustics  study  for  the 
Bermuda-Eleuthera  range  is  obvious. 

Future  analyses  using  the  BEAR  Buoy  data  set  should 
examine  in  detail  the  modal  structure  and  changes  in  dis- 
tribution by  mode  over  the  period  or  record. 
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